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Abstract. The definite integral 



M(a) 



r*/ 2 x 2 dx 

x 2 + ln 2 (2e~ a cosx) ' 
is related to the Laplace transform of the digamma function 



4 r 

7T Jo 



e- aa ip(s + 1) ds, 





by M(a) = L(a) + ■y/a when a > ln2. We establish an analytic expression for 
M(a) in the complementary range < a < In 2. 



1. Introduction 

The classical table of integrals by I. S. Gradshteyn and I. M. Ryzhik [7] contains 
a large collection organized in sections according to the form of the integrand. 
In each section one finds significant variation on the complexity of the integrals. 
For example, section 4.33 — 4.34, with the title Combinations of logarithms and 
exponentials, presents the elementary formula 4.331.1: for a > 0, 

7 + In a 



(1.1) 

where 7 is the Euler constant 



n 1 

(1.2) 7= lim > Inn, 



n — >oo £ — ' h 

fe=l 



as well as the more elaborate 4.332.1 and 4.325.6 

\nxdx Z -1 , , fl\ dx 2tt /5 , , (\ 

= / In n - = - n 2tt - n V - 



, In In - -= = -= -bi27r -lnr - 

e~ x -l J \x)x 2 -x + l V3 V 6 V 6 



The difficult involved in the evaluation of a definite integral is hard to measure from 
the complexity of the integrand. For instance, the evaluation of Vardi's integral, 

(1.3) f"\^n xdx= [\ b (l)j^ = lJ T ^\ 



tt/4 



1 + ^2 2 I r(i) 



that appears as 4.229.7 in [7], requires a reasonable amount of Number Theory. The 
second form is 4.325.4, found in the section entitled Combinations of logarithmic 
functions of more complicated arguments and powers. The reader will find in [15] 
a discussion of this formula. 
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It is a remarkable fact that combinations of elementary functions in the integrand 
often exhibits definite integrals whose evaluations are far from elementary. We have 
initiated a systematic study of the formulas in [7] in the series [T} [21 [HI EH [HI Q~2] • 
The papers are organized according to the combinations appearing in the integrand. 
Even the elementary cases, such as the combination of logarithms and rational 
function discussed in [2] entail interesting results. The evaluations 

<"> jfuSpr = i[i-(i + »)-]M»4Mi + !» 



1 1 / 1 

1 ^ 1 / n — 1 



n (l + j\ \ j 



\s(j + 1,2)\V, 



for b > and n G N produces an explicit formula for the case where the rational 
function has a single pole. Here s(n, k) are the Stirling numbers of the first kind 
counting the number of permutations of n letters having exactly k cycles. The case 
of a purely imaginary pole is expressed in terms of the rational function 

(1.5) Pn(z)=£- 



as 



with 



in tdt _ ( 2 :) 



(l + f 2 )^ 1 2 2 » 



^-4 tan 1 x+p k (x) 

go(x) +pn{xj hi 



^ 2fc + 1 



(1.6) go (x) — hi £ tan 1 x — / dt. 



(1-7) , m = -2~ 2 "' 




The special case x = 1 becomes 
Z" 1 In i eft 

y (i+i 2 ) n+i 

where 
(1.8) 

is the Catalan's constant. The values 

do not admit a closed- form (in the sense of [T3]), but they do satisfy the three term 
recurrence 

(1.10) (2k + 1W(1) - (3k + l)p fc (l) + fep fe _i(l) = 0. 

The study of definite integrals where the integrand is a combination of pow- 
ers, logarithms and trigonometric functions was initiated by Euler [5], with the 
evaluation of 

«r/2 n 

(1.11) / x \n(2cosx)dx = C(3), 

Vo 16 
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and 

r-7r/2 



C(3), 



(1.12) / x 2 ln(2 cosx) <ir = 

that appear in his study of the Riemann zeta function at the odd integers. These 
type of integrals have been investigated in [8] , [16) . The intriguing integral [3] , 



(1.13) 



/ a; 2 In (2 cos x) dx = (U) 

Jo 16 



1440' 



was first conjectured on the basis of a numerical computation by Enrico Au-Yueng, 
while an undergraduate student at the University of Waterloo. 

Recently O. Oloa [13] considered the integral 

^ 2 x 1 dx 



(1.14) 



4 f 
M(a) := - \ 
it Jo 



x 2 + In (2e a cos a;) 



and later established the value 

4 rl 2 x 2 dx 1 
(1.15) M(0) = - / , , . = - (1 + ln(27r) - 7 ) 



7T 



x 2 + In (2 cos a;) 2 



Oloa's method of proof relies on the expansion 

(1.16) % = xsm2x+Y{-\) n - x ( ^ - -^±l- )xsm(2nx) ) 

1 7 x 2 + ln 2 (2cosx) ^ ' \n\ {n + l)\ ) V ; ' 

where 

(1.17) a n := f (t) n dt, 

Jo 

with (t) n = t(t + 1) • • • (t + n — 1) the Pochhammer symbol. The standard relation 

(1.18) (t)„=]T| S (n,fc)|i fc , 

fe=l 

gives 

|s(n,fc)| 



(1.19) = £ k + 1 • 

fc=i 

M. L. Glasser and D. Manna 6] introduced the function 

(1.20) L(a) := / e" QS V(s + 1) da, 

Jo 

where ^(a;) = 4=- lnT(x) is the digamma function. Integration by parts and using 
(jl.ip gives 

(1.21) L(a) = -7-lna + a / e _a * lnT(t) dt. 

Jo 

The main result in [6] gives a relation between M(a) and L(a). 
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Theorem 1.1. Assume that a > In 2. Then 

M{a) = L(a) + -. 

That is, for a > In 2 , 

/>OC 

(1.22) M(a) = - -7-lna + a / e _a * lnr(i) <2t. 

a Jo 

The proof in [B] begins with the representation 



,tt/2 

(1.23) / cos" x cos ax dx — 

Jo 



ttI> + 2) 



2»+ 1 (v + 1) r(i + f + f ) r(i + f -§)' 



borrowed from 3.621.9 in [7]. Differentiating with respect to a, evaluating at a = s, 
and using "0(1) = — 7 yields 

2S+2 rrc/2 

(1.24) ?/;(s+l) = / xcos s x sin(sx) dx — 7. 

I" Jo 

Replacing in (|1.20|) produces 



7 4 



a 7r 



00 



"'O 



(1.25) L(a) + - = --Im / / xe s(ln [ 2e_a cosx h™) dxds. 



The identity (1 1 . 22[) now is an immediate consequence of the s-integral: 

(1.26) / e ^H 2e ' acosx }-"ds= = r. 

J ix — In [2e~ a cos x\ 

The authors also succeed in a series expansion of M(a) while they recognize as 
a hypergeometric function in two variables, and noted (quoting from [6]) strongly 
suggests that for general value of a, no further progress is possible. This hypergeo- 
metric interpretation led the authors [6] to 

1 r 1 

(1.27) M(0) - 1 + - / i(l - t) 3 F 2 (1, 1, 2 - t; 2, 3; 1) df 
for which they invoke 

(1.28) 3 ^l.l,2-,;2,3;l)= 2 "-\7f + 1 >> 

to enable them demonstrate a new proof of (|1.15|) . 

The graph of M(a) shown in Figure 1, obtained by the numerical integration 
of (|1.14|) . has a well-defined cusp at a = In 2. In this paper, we provide analytic 
expressions for both branches of M(a). The region a > In 2, determined in [B], is 
reviewed in this section. The corresponding expressions for < a < In 2 will be the 
content of the next section. 

2. The case < a <ln2 
Our starting point is the identity 

(2.1) M(a) = -— Im / e~ at f ^±^1 dxdt 
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2.0 r 



1.5 - 



0.0 — 
0.0 



FIGURE 1. The graph of M(a) for < a < 2 



as illustrated in [5J. We outline the proof here for sake of the reader's convenience. 
The identity 

T. T 

(2.2) Im ■ 



yields 
(2.3) 



ix + In [2e~ a cos x] x 2 + In 2 [2e- a cosx] 

4 r 

Mia) = - Im 

n Jo 



; dx 



ix + In [2e~ a cos x] 



Under the assumption a > In 2, we have 



(2.4) 

which implies 
(2.5) 



3 s ln[2e a cos cc] -\-ix 



ds 



In [2e a cos x] ' 



7T/ Z /> OC 



r/2 

M(a) = — Im / / .)•( 

T J -tv /2 Jo 



s(ln^2e a cosa;j+'tx) 



da; ds, 



where one uses the fact that the imaginary part of the integrand is an even function 
of x. One more identity 

^2 6) e isx . e sln [ 2e " cosx ] — e s i n [ e_ °( 1 + e2 ")] 

and the change of variables x i— > x/2, give reason to 

(2.7) M{a) = —lmj f ie sln l r ' (1+e ")] dsdx. 

Then evaluate the s-integral to obtain 

1 , r xdx 



(2.8) 
The formula 
(2.9) 

now gives QUI]) from |[2T5)) . 



27r Im ,/_ 7r ln[e- Q (l + e")]' 



1 f 1 u*dt 
In u 7 u — 1 
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Note 2.1. Even though the proof outlined here is valid for a > In 2, the identity 
pOp holds for a > 0. 

Notation: we use b = e° — 1 and assume < a < In 2, so that < b < 1. 
Expanding the terms (1 + e™)' and 1/(1 — be~ lx ) in power series produces 

-1 ^ oo oo 



M(a) = -\ 



- / / xe- at VVi 5 sin[a;(fc - j - 1)] dx dt. 

77 JO J -it ,-_nt_n V*/ 



,j=0 fc=0 

The term corresponding to fc = j + 1 disappears and computing the x-integral we 
arrive at 

(2.10) M(a) = e Q / e-"^^(-lp fc | rfi 



j + 1 - k 

j=0 k=0 J 



oo 



3 

Lemma 2.1. Let t e K and j e N U {0}. Then 

(-i)V t \ 

Proof. The integral representation (3.268.2 in [7]): 



(2.H) ^^^^.j = g Wi+ i)_^ + i)]. 



(2.12) ^(p+i)_^( g + i) = _ / 

Jo 1 — x 

yields 



1 X p - xi 



(2.13) iP(p+l)-i>(q + l)=J2(- i y~ 1 



P\ 1 



Therefore the result is a consequence of the identity 

'tY 1 y, (-1)"' / t \ y, / t \ _ * (-l)m-l ^ /, 

^— ' m \m + k J ' m ' 

m— 1 v 7 m—1 v 7 m—1 

Apply the difference operator Aa(fc) = a(fc + 1) — a(fc) and use 

< V 1 / i A ft\~* f t \ m ( t \~ 1 / i+1 



v fc + ly \m + /c + ly y/cy \m + k J k + l\k + lj \m + k + 1 
to write the derived equation as 

(2.15) _^iLfv_i r ( ' + 1 LaVH 

m—1 v 7 m—1 v 

The left hand side of (|2.14|) reduces to — l/(fc + 1) in view of the classical identity 

t+1 \ _ / £ 
v m + fc + 1/ ~~ \A; + 1, 



(2.16) E^ 1 )™" 1 



A simple evaluation of the right hand side in (|2. 14[1 also produces — l/(fc + 1). We 
conclude that, up to a constant term with respect to the index fc, both sides of 
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()2.14|) are equal to the harmonic number Hk- The special case k — shows that 
this constant vanishes. □ 

Continuing from (|2.10p , we thus have 
(2.17) 

j=0 k=0 



M { a) = e"[ E W E TTlTfc + T [ '"T^W 1 )* 
Jo j=0 fc=0 J + 1 - * o Jo j=1 \J/ 

a pi 00 / f \ 

r/ .-E"G)* (,+1)A 



To simplify the first term in the previous expression observe 

hl (l + & ) - /A ae q« 
b ^\k) b ■ 

We deduce that in ([2~T7| the first term is a/(l - e~ a ). 

The reduction of the second term in (|2.17p employs the following result: 
Lemma 2.2. Let < a < In 2 and tel. Tften 

(2.18) ^ e- Qt E QJV'C? + 1) * = ln(l - O + ^ — dt. 

Proof. The Stirling numbers s(j, fc) satisfy 

so that 



(2-19) j!(f) = $>(?, 



fe=0 



(2.20) / e- ai E^( .W0' + 1)#=- -+e- a E(^ +lc ^-^-iMi + 1 )> 

•/O „_n V/ a „-_i 



with 



j=0 j=l 



( 2 - 21 ) a i ( a ) = ~j ~FTT 



1! ' — ' a" 

J k=o 



The result now follows by summation by parts and the identity 

\s(j,k)\V ln fe (l + 6) 



(2.22) E 



jl k\ 



Therefore, the second term in (|2. is 

lnfl - e~ a ) 1 f 00 e~ at 
(2.23) second term = — H / dt. 

K 1 1 - e- a 1 - e- a J 1 t 



□ 
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Finally, the third term in (|2.17j) is 



, • , e 
third term = — — 
o 



A direct computation shows that J Q ^>(f + 1) dt — 0, and integration by parts 
gives 

(2.24) third term = - _ a / e _a ' lnT(< + 1) dt. 

The identity lnT(< + 1) = lnT(t) + hit now yields 

third term = - °_ a ^ e" Qt ln<dt + jf e" at lnT(i) . 

Replacing (2.17), (2.19) and (2.21) into (|2.17|) provides the following expression for 
M{a): 

a 7 ln(l - e~ a ) 



(2.25) M{a) 



l-e~ a a 1 - e~ a 



/ e- at lntdt + / e- at \nT(t)dt. 



The term j/a comes from the index j = in the sum (|2.18p . Next, we make use 
of (jl.ip to state our main result which, incidentally, is complementary to Theorem 
11.11 This settles a conjecture of O. Oloa stated in [T3] . 

Theorem 2.1. Assume < a < In 2. Then 

, . 7 a + ln(l — e~ a ) — 7 — lna a f 1 . 

Ma = 2 + =— ; ^ + 1 / e- at In r t 

a 1 — e° 1 — e a 



This result enjoys a form similar to Theorem ll.il 
Corollary 2.1. Assume < a < In 2. Then 

MM ^ + ( a+M1 -e : ) + r(o, a )) + _j_| e _. V(t + 1)a 

where r(0,a) is i/ie incomplete gamma function. 

Proof. Split up the first integral in ()2.25|) and integrate by parts. □ 
Differentiating (|2.1|) with respect to a at a = 0, and use the classical value 

(2.26) / \nT(t)dt = -ln27r 

Jo 2 

and also 

(2.27) / tlnr(t)dt = + iln2?r- — 
1 ' Jo 2tt 2 6 12 

obtained in [4], produces 

r /2 x 2 ln(2cosx)dx 7tt tt ln27r C'(2) 
l2 ' 2M Jo (x 2 + In 2 (2 cos x)) 2 ~ 192 + "96 167' 

Further differentiation of (|2.1[) produces the evaluation of a family of integrals 
similar to (f2~2"g)) . 
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The integral in (|2.1| can be expressed in an alternative form. Define 



( 2 - 29 ) A w=i Era- 111 



Observe that A(0) = 7, so A(z) is a generalization of Euler's constant. 
Lemma 2.3. Let a > 0, c = 1 — e~ a and define A = ln27r + 7. TTiera 



(, 30 ) /'e-mrw,^^)-^^)^^ 



_ x a 2 + 47r 2 j 2 ' 



Proof. Expand the exponential into a MacLaurin series and use the values 

'n + l\ (2ft)! 
v 2fc- l) k(2n) 2k 



L-J 

- ^B-)*(" 2 r)^^ +1 ) 



71+ 1 



given as (6.14) in 0j. Then calculate by interchanging the resulting double sums. 

□ 

The next corollary follows from the identity M (a) = L(a) + 2. 

Corollary 2.2. Assume < a < In 2 and Zei c = 1 — e~ a . Then 

f°° „ f , , x , 7 + In a yl(a — c) 1 , / a \ In j 
2.31 / e- Qt lnr (i = -- + ^ — - — —A — +2V — j-—. 

Jo ace a a 2 c 2a \2irJ a 2 + An 2 ] 2 

We now state two identities involving the function f(t) = 2 _t lnT(i). The proof 
of these identities was supplied to the authors by O. Espinosa. 

Lemma 2.4. The identities 

(2.32) l°°f(t)dt = 2^ 1 /(^- 7+ 1 ^ 2 ln2 

f°° „ , Z -1 , . , . (7 + lnln2)(l + 21n2) -1 

(2.33) / */(t)dt = 2/ t + i/p)*-!! Ji^ > 

Jo Jo In 2 

Proo/. The function f(t) satisfies f(t + l) = \f(t) + \2~ l Int. Splitting the integral 

(2.34) / f(t)dt= / /(t)dt + / f(t + l)dt 
Jo Jo Jo 

and using (jl.ip gives the first result. The proof of (|2.33|) is similar, it only requires 
differentiating (jl.ip with respect to the parameter a. □ 

The reader will check that (|2.32[) is equivalent to the continuity of M(a) at 
a = In 2. The identity (|2.33[1 provides a proof of the last result is worthy of singular 
(pun-intended) interest. 

Theorem 2.2. The jump of M'(a) at a = In 2 is 4. 
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The integral 



4 I 
M(a) :=- - 
n Jo x 



3. Conclusions 

4 r /2 x 2 dx 



satisfies 



/o x 2 +ln 2 (2e a cosx) 

(3.1) M(a) = -+ [ e- at ip(t + 1) dt, 

a Jo 



for a > In 2 and 



. . 7 (a + ln(l - e- a ) +T(0,a)) 1 f 1 „, ,. 

M(a) = 2 + i f — A + — ^ / + 1) dt 



for < a < In 2. 
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